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Abstract 

We investigate the leptogenesis with almost degenerate neutrinos, in the frame- 
work of democratic mass matrix, which naturally explains the large mixing angles 
for neutrino oscillations as well as quark masses and mixing matrix. We find that 
the baryon asymmetry in the present universe is explained via the decays of right- 
handed neutrinos produced nonthermally by the inflaton decay. The model predicts 
neutrinoless double beta decays accessible in near future experiments. 



One of the fundamental problems in particle physics is to understand the observed 
masses and mixing angles of quarks and leptons. There have been, in fact, proposed 
many models of the mass matrices for quarks and leptons. Among them the hypothesis 
of democratic mass matrix [|IJ is very interesting not only theoretically, but also phe- 
nomenologically. Here, the three families of quarks and leptons are treated in an equal 
footing, i.e. a permutation symmetry S3XS3 is imposed among three families. The demo- 
cratic mass matrix has been known very successful in explaining the observed masses and 
mixing angles for quarks 0. This approach has been extended to the charged lepton 
and neutrino sectors, and it has been found that the observed large mixings of neutrinos 
can be rather easily understood due to the nearly diagonal neutrino mass matrix [§, |J.[] 
Surprisingly, the model also accounts for the small mixing angle U e 3, required from the 
CHOOZ experiment ||, although other two mixing angles are large [[5], §. Furthermore, 
this scheme yields almost degenerate Major ana masses for neutrinos ant it predicts neu- 
trinoless double beta decays accessible in near future experiments, such as GENIUS |J, 
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In this letter, we discuss the leptogenesis [|1J] in the framework of the democratic mass 
matrix. We show that the baryon asymmetry in the present universe can be generated via 
the decay of heavy right-handed neutrinos. Here, the right-handed neutrinos are assumed 
to be produced nonthermally through inflaton decay |15 . 

Democratic mass matrix 

Let us start by giving a brief review on the work in Ref. ||, which imposes a permu- 
tation symmetry Ss(L)x S3(i?) on the three families of left-handed lepton doublets (lu) 
and right-handed charged leptons (em). We assume that In transform as 3l = 2l + 1l 
under S 3 (L), and e Ri transform as 3# = 2 R + 1 R under S 3 (i?). Then the mass matrix 
for charged leptons which is invariant under the S 3 (L)x S 3 (-R) is uniquely determined as 
follows: 



M, 



(o) 



Mi 
3 



1 1 1 
111 
1 1 1 



(1) 



This mass matrix is often called "democratic" mass matrix [|TJ. Two mass eigenvalues of 
this mass matrix vanish. Thus, we introduce a small symmetry-breaking mass term as 



1 Rcf. Q considers Dirac or Majorana neutrinos. Thus, their model does not have a group theoretical 
reason for the family-diagonal and degenerate neutrino masses. 
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perturbations with a diagonal mass matrix 



M, 



(i) 



-ei 

El 

o o <y, 



(2) 



Then the matrix Mi = + is diagonalized as 

V/M^diagK, 774,774) 



where 



7771 



Mi 



'Si 6 



777 Q 



777o 



Mi 1 + 



with 



6 = 2 (a? + 3^ 



1/2 



(3) 



(4) 



(5) 



where terms with higher order of 5i and ei have been ignored. The unitary matrix Vi 
which diagonalizes the mass matrix M\ is given by VJ = ABi, where 



and 



B, 



A 



cos d\ — sin Qi 
sin 6i cos 
A/ sin 39i \i cos 



1/V2 1/V6 1/V3 
-1/V2 1/^6 1/V3 
-2/V6 1/V3 



(6) 



-A/ sin 28 1 
-A/ cos 2^ 
1 



tan 28, 



y/3Ei 



A, 



3y/2' 



(7) 



m e , m: 



and rng 



777 T J are 



The empirical masses for the charged leptons {m\ 
given by adjusting appropriately the parameters Mi, Ei and Si as 

Mi = 1719 [MeV], e x = 0.00665, 5i = 0.095 , (8) 

which leads to sin^ ~ 0.066. 

Let us turn to the neutrino sector. Before introducing heavy right-handed neutri- 
nos, we discuss the mass matrix for the light Majorana neutrinos in the effective non- 
renormalizable operator: 



fij 

2M 



(9) 
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where H denotes the Higgs doublet. There are two independent mass matrices that are 
invariants under Ss(L) symmetry ||: 



(10) 



Here we take the first form for the time being, since it automatically leads to the large 
mixing angles for both of the solar and the atmospheric neutrinos. (We will add the 
second mass matrix later.) Then, the neutrinos are degenerate in mass. With symmetry 
breaking terms in diagonal elements, the neutrino mass matrix is given by 
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The neutrino-mixing matrix in the basis where the mass matrix for the charged leptons 
is diagonal is given by U = (Vi)' as R 



U = (AB,) f ^ A ] 



l/y/2 -l/y/2 
1/VE 1/y/E -2/VQ 
l/y/3 l/y/3 l/y/3 



(12) 



As denoted, this indicates nearly bi-maximal neutrino oscillations, i.e. 

8 



sin 2 29 



12 



sin 2 26 



- ' ~ 9 ' 



(13) 



where 9\2 and #23 are the solar and the atmospheric neutrino mixing angles, respectively. 
Furthermore, the U e 3 element automatically vanishes, U C 3 — 0. The mass squared dif- 

2M 2 u e u and 



ferences for solar and atmospheric neutrino oscillations are given by Am. 

2M. 2 U 5 V1 and hence the breaking parameters should satisfy e v j 8 V 



solar 



Am atm - 



Am sola, 



/Am 



2 

atm" 



The experimentally observed values for these quantities are given by 0, [16j 



llul ^{(1.0-6.0), 3.2}xlO- 3 eV 2 
4i ar = {(2-50), 4.9}xlO~ 5 eV 2 



Am 2 
Am 

sin 2 2^3 
tan 2 ^i2 



= {(0.83-1), 1} 
{(0.2-0.8), 0.37} 



(14) 



Here, we have neglected the small mixing angle 81 of charged leptons. 
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where the last numbers in the parentheses denote the best fit values. Here, we have taken 
the large mixing MSW (LMA) solution to the solar neutrino problem, which is the most 
preferable in the present experiments ||. 

Recently, a positive indication for neutrinoless double beta decay has been reported, 
with v e -v e component of the neutrino mass matrix |j7i„ e „ e | ~ 0.05-0.84eV |17|. Therefore, 
it might be interesting to consider the degenerate neutrinos of mass of order O(0. leV). 
Hereafter, we take ~ M. v ~ O.leV as a representation. In this case the breaking 
parameters in the neutrino mass matrix are obtained as e u = Am s 2 olar / (2.M 2 ,) ~ 0.001- 
0.025 and b v = Am^ tm /(2.A/^) — 0.05-0.3. We see that their orders of magnitudes are 
just around those of the breaking parameters in charged lepton mass matrix, Ei = 0.00665 
and 5i = 0.095, in Eq. (§). 

One may worry about that the predicted value of the mixing angle for the solar 
neutrino #12 in Eq. ( |13|) is too large to fit the experimental value. Even if we include the 
effect of nonzero sin#/ ~ 0.066, the mixing angle of the solar neutrino oscillation is only 
slightly reduced to sin 2 2# 12 — 0.99. However, a deviation from the nearly maximal mixing 
for solar neutrino oscillation can be easily implemented by introducing small off-diagonal 
elements k v in the neutrino mass matrix Eq (|TTJ) : 



M v = M, 



1 K v K v 

K u 1 + E v K v 

l + 5 v 



(15) 



Notice that the matrix of the off-diagonal elements has the form of the second matrix in 
Eq. ([TO]), which is also invariant under the Ss(L) as noted before. The mixing matrix of 
neutrino-oscillations is then given by U = (Vi)*V v . Here, V v is defined as 

V}M V V V = diag^, m w 2 , m») , (16) 

where ~ m v 2 ~ M. v (l + e u ), m\ ~ M. V {X + 8 U ). Here and hereafter, we require 
k v < e u in order to ensure the nearly bi-large mixing for neutrino oscillations. 

The neutrino mass matrix given in Eq. ( |T5D has another interesting prediction on the 
size of the U e 3. When the neutrino mass matrix is diagonal as in Eq. (|II|), the neutrino 
mixing matrix U is entirely determined by Vi given in Eqs. (](]) and ([I]). In this case, we 
can immediately obtain the U e 3 element as 

2 A; 

U e3 = -= sin 9i + -j= sin 30j ~ -0.05 . (17) 

V6 v3 
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In the case of the general neutrino mass matrix, the above prediction on the U e 3 is shifted 
by the non-zero elements of the V v . The leading deviation of the U e 3 from the value 
presented in Eq. (fi7|) is given by 

cos 9[ 



AU 



c3 



V2 



z/13 



V,. 



u23) 



(18) 



We now consider the specific case where the neutrino mass matrix is given by Eq. flT5| ) 
One might naively imagine that the deviation of the U e 3 is given by 

' K 



AU 



<3 



o 



S v 



10 -d-2) 



(19) 



from the above general argument. If this is true, we can have no precise prediction on 
the size of the U e 3. Interestingly, however, this is not the case. Because of the symmetric 
form of the neutrino mass matrix in Eq. (|l5l), the leading contribution given in Eq. (|18|) 
beautifully cancel out. Actually, the leading terms in V v \z and V^z are given by 



Kl3 

K23 



1± 

K,. 



+ 



v v 



Su si 



(20) 



Hence, the deviation of the U e 3 in this model is highly suppressed, AU e 3 ^ 10~ 3 , which 
allows us to have a precise prediction, |U e3 1 ~ 0.05.0 This will be clearly shown in the 
numerical calculations presented in this letter. In the future, this prediction on the size 



of the U e3 will be tested in the long baseline experiments, such as JHF [18 



The most natural way to explain small neutrino masses is the "see-saw" mecha- 
nism ULy[|. In this case, the mass matrix M v for the light neutrinos is given by 



My 



m. 



vD M R m »D 



(21) 



after integrating out the heavy right-handed Majorana neutrinos Nm. Here, Mr and m V D 
are the mass matrices of the right-handed Majorana neutrinos Nja and the Dirac mass 
for In and N^j, respectively. In terms of the Yukawa couplings, the Dirac neutrino mass 
matrix is given by {m V D)ij = (hr>)ij(H), where (ho)ij is defined as C = (hr^ijNmlLjH. 
The S3 invariant matrices of Mr and hp are uniquely determined as 



Mr = Ail 



1 K R K R 
KR K-R 1 



h 



D 



k 



D 



1 K D K D 

Kb 1 k d 



(22) 



In fact, V„i3 = V V 23 in the limit of the vanishing e v . 
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where we have assumed that the Nm transform as 3^ = 2^ + I/, under the Ss(L) symme- 
try.0 One can easily imagine that the assumption of kr ~ ko(~ k v ) for the off-diagonal 
elements is the simplest and natural way to obtain the mass matrix for the light neutrinos 
in Eq ( pT5|) with e u = 5 U = 0. In order to induce the symmetry breaking parameters e v 
and 8 V in the neutrino mass matrix in Eq. (|15|), we introduce breaking parameters Ed and 
5d in diagonal elements of Yukawa matrix ho, as in charged leptons: 



h D 



1 D 



1 K D K D 
K D 1 + E D K D 
K D K D 1 + 8 D 



(23) 



Then, we can obtain the required mass matrix of the light neutrinos by taking the following 
pattern of the perturbations: 

k r ~k d (~k v ) < e d (~ e„) ~ 0(0.01) < 8 D (~ <5„) ~ O(O.l) . (24) 

The above breakings of the S3 symmetries might be understood as follows. Suppose 
that the S3 breakings are originated only from the Yukawa coupling of Higgs field H, as in 
mass matrices for quarks and charged leptons. Then, the Ss(L) symmetry is broken in the 
neutrino Dirac Yukawa matrix by the parameters Ed and 5d as in charged lepton sector, 
since it comes from the coupling of Higgs field. On the other hand, the mass matrix of the 
right-handed neutrinos M R in Eq. (p2|) is assumed to be S 3 invariant since it is decoupled 
from the Higgs field H. 

The mass matrix Mr for the heavy right-handed neutrino leads to two exactly de- 
generate right-handed neutrinos and a slightly heavier /lighter one, with masses M Ri = 
{A4r(1 — Kr), A4r(1 — Kr),A4r(1 + 2k,r)}. As we will see in the remainder of this letter, 
the off-diagonal elements kr and kjj, which are required to explain the neutrino oscillation 
experiments, play a crucial role in the leptogenesis. 



4 If Nm transform as 3^ = 2r + 1^ under the Sa(i?), the Dirac Yukawa matrix becomes democratic 
one as in Eq. (Q), which invalidates the almost diagonal mass matrix for light neutrinos in Eq. (|l5|). 
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leptogenesis |l4j 

Now let us discuss the leptogenesis in the present model. Notice that the conventional 
leptogenesis scenario where the right-handed neutrinos are produced by thermal scatter- 
ings after the inflation is somewhat difficult in the case of the degenerate neutrinos of 
m\ — m 2 — m 3 — C(0.1)eV. This is because the out-of-equilibrium condition cannot 
be satisfied in this case and the amount of produced lepton asymmetry is strongly sup- 
pressed pOfl . (See Appendix.) Therefore, we consider the leptogenesis via decays of the 
right-handed neutrinos Nm which are produced non-thermally. A natural mechanism of 
such a non-thermal production of Nm is the decays of the inflaton ip into the Nm . 
Hereafter, we will consider the supersymmetry (SUSY) theory, although the following 
discussion does not change much in the non-supersymmetric case. 

If CP is not conserved in the Yukawa matrix ho, the interference between decay 
amplitudes of tree and one- loop diagrams results in the lepton-number production [Q. 
The lepton-number asymmetry per a decay of right-handed neutrino N Ri is given by |T3L 



Sj r(iv ffi -> i Lj + h)- Ej r{N Ri -> H~ + ff) 
Ei r(iv ffi -> i Lj + h) + Ej r(N m -^TTj + H) 



87T (h D h ] D 



J2im[{(h D hi) ik y 



M} 



'Ml 
M? 



(25) 



I'D) H k^i 

where N Ri , l^j, and H (l^ and H) symbolically denote fermionic or scalar components of 
corresponding supermultiplets (and their anti-particles), and Ty{x) and J-s{x) represent 
the contributions from vertex and self-energy diagrams, respectively. In the case of the 
SUSY theory, they are given by [p2j] 



x 



(26) 



Here, we have assumed that the mass difference of the right-handed neutrinos is large 
enough compared with their decay widths, so that the perturbative calculation is ensured. 
(We will justify this assumption later.) 

In the present model, the masses M Ri are almost degenerate, M Ri = {M. R {1 — 

K R ), A4 R (1 — Kr),A4r(1 + 2k r )}. Thus, x ~ 1 and the self-energy contribution Ts(x) is 

much larger than the vertex contribution jFy(x).Q In the leading order in perturbation, 

5 If we take into account the effect of the finite decay width, Ts{x) vanishes for a; — > 1. Thus, 
T s {Ml/Mf) in Eq. (||) vanish for % = 1, A; = 2 and i = 2, k = 1. See also footnote §. 
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the asymmetry parameters are given by 

12tt D Re(K R ) D ' 9vr D Re(n R ) D ' 6 K J 

where we have neglected higher order terms in the expansions of Ed, 5d, Kq, kr and 
Ed/Sd- As for the CP phase, we have assumed that the complex phase exists only in the 
off-diagonal elements and kr for simplicity, and have taken the other parameters to 
be real. Notice that all the decays of Nm generate the lepton asymmetry with the same 
sign, namely, they contribute in a constructive wayfl 

The ratio of the lepton number density til to the entropy density s produced by the 
inflaton decay is given by [|TjJ 

- = !£^-, (28) 

where Tr is the reheating temperature after the inflation, the mass of the infla- 
ton, and the branching ratio of the decay channel of the inflaton to N Ri , i.e., 
= B r ((f) NffiNjn). Here, we have assumed that the inflaton decays into a pair of 
right-handed neutrinos, and Mm > T# in order to make the generated lepton asymmetry 
not washed out by lepton-number violating processes after the iV^'s decay. Notice that 
the inflaton mass should satisfy > 2M. R in order to make the decay ((f) — > N^Nm) 
kinematically allowed.[] After being produced, a part of the lepton asymmetry is immedi- 
ately converted |14]] into the baryon asymmetry via the "sphaleron" effect ||24|| , since the 
decays of N Ri take place much before the electroweak phase transition: 

= , (29) 
s s 



where C is given by C ~ —0.35 in the minimal SUSY standard model (MSSM) p5| . 
Therefore, the amount of the baryon asymmetry in the present model is estimated as 

M ^ 0.3 x 10- 10 x (£<?) +B?>) ( Tr ^ (2MR ^ ( M » ^ f6D ^ lm(KD 



V10 8 GeV/ V m 4> J VO.leV/ V - 1 / Re(«B) ' 
(30) 

6 One might wonder if this fact might conflict with the argument that the generated lepton asymmetry 
must vanish in the limit of exactly degenerate masses. However, if the mass differences of the right-handed 
neutrinos become smaller than the decay widths of them, the perturbative formula !Fs(x) in Eq. (|2^) no 
longer holds and we should take into account the effect of finite widths of Njh. Actually, it was shown 
that Ts{x) vanishes in the limit of exactly degenerate right-handed neutrino masses if the finite decay 
widths are appropriately taken into account |p3[ . 

7 In this letter, we assume a perturbative decay of the inflaton. 
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where we have used the see-saw formula M. v ~ k 2 D (H) 2 / M.r^ We see that the empirical 
baryon asymmetry ~ (0.4-1) x 10~ 10 |26| is obtained with a reheating temperature 



Tr ~ 10 8 -10 10 GeV, for a natural choice of the parameters, say, + ~ 0.1-1, 
2M R ~ (O.l-l)m^, 5 D ~ 0.1 and Im(/t D )/ Re(K R ) ~ 1. 

Notice that the mass differences of the right-handed neutrinos, AMr ~ 3Ke(K,R)AiR 
is much larger than the decay widths of them, — (k^/ Air) Mr, as long as kr ^> 

fe|,/127r ~ 10~ 6 x (Mv/0.leV)(M R /l0 10 GeV), so that the perturbative formula F s {x) in 
Eq. (|26|) is justified. The assumption > Tr is also easily satisfied for M& > 10 10 GeV. 

In Fig. m, we show the histograms of the ur/s, [U e 3 1 and for the LMA solution. 
Here, we randomly generate the small perturbations as in Eq. (p4|) and collect the data 
set of Ub/s, I U e 3 1 and M. v if the generated mass spectrum and mixing angles satisfy the 
conditions obtained from the neutrino oscillation experiments given in Eq. (0). Then 
we plot the frequencies of these quantities in the vertical axes. We also present a plot 
of the obtained mass differences and mixing angles for the solar neutrino oscillation in 
{tan 2 6*i2, Am s 2 olar } plane. We take the range of the perturbations as follows: 

5 D = (0.25 - 1.5) x 10" 1 , e D = 5 D 10~ (a3 - 2 - 5) , 

KD = £ D 10-(°- a V(°- 2 ^, kr = \K D \e«-*~fr, (31) 

where we take Ed and 5d to be real, for simplicity. If we randomly generate the perturba- 
tions within the above range, the neutrino mass spectrum and mixing angles satisfy the 
conditions given in Eq. (O) with the rate of nearly 20%. The baryon asymmetry ur/s is 
calculated with the normalization 

One can see that the analytic estimation explains the numerical result quite well and that 
the required amount of the baryon asymmetry can be easily obtained in the democratic 
model with the natural scale of perturbations. Another interesting prediction can be seen 
from the upper-right figure. The amplitude of U e3 is accurately predicted as |U e3 | ~ 0.049, 
as denoted in the first part of this letter. 

In Fig. we show the histogram of the the ur/s, [U e3 1 and Ai u for the LOW solution. 

8 In the MSSM, (H) = 174GeV x sin/3, where tan/3 = (H)/{H') and H' is the Higgs field which 
couples to the down-type quarks (and charged leptons). In Eq. (|30|), we have taken sin/3 ~ 1. 



9 



In this case, we take the range of perturbations as follows: 

S D = (0.25 - 1.5) x HT 1 , e D = S D 10^ 37 ~ 5 - 8) , 

KD = e D 10-<f>-°'%W-*>, k r = \K D \e^^\ (33) 

Other conventions are the same as those in Fig [l]. The resultant baryon asymmetry is 
almost the same as in the LMA solution, which is easily understood from the Eq. (0). In 
the case of the LOW solution, the off diagonal elements in the light neutrino mass matrix 
are much smaller than those in the LMA solution k v ~ e u ~ 5 U x 10~^~ 6 \ which results 
in the much more precise prediction on the amplitude of the U e 3. 

Conclusions 

In this letter we have investigated leptogenesis with almost degenerate neutrinos, in 
the framework of democratic mass matrix, which explains very successfully the observed 
large mixings of the neutrinos as well as quark masses and mixings. The almost degenerate 
Majorana neutrinos with masses of order ~ C(0.1)eV induce a considerable rate of 
the neutrinoless double beta decays, which is accessible in near future experiments. We 
have shown that the empirical baryon asymmetry is well explained by the decays of the 
right-handed neutrinos produced in the inflaton decay. 

In this model, the U e 3 component of the mixing matrix for neutrino oscillations is 
predicted as | U e 3 1 — 0.05, which is a direct consequence of the nearly S3 symmetric form 
of the neutrino mass matrix. Such a value of U e 3 will be also tested in the long baseline 
experiments, such as JHF [TH . 



Appendix 

In this appendix, we show that if the light neutrinos are degenerate in masses of 
~ O(0.1)eV, the produced lepton asymmetry is highly suppressed in the case of the 
leptogenesis scenario where the right-handed neutrinos are produced by thermal scatter- 
ings. We should stress that the following discussion is a generic consequence of the almost 
degenerate light neutrinos, independent of the specific models. 

The wash-out effect of the lepton asymmetry in the leptogenesis by the decays of ther- 
mally produced right-handed neutrino Nm crucially depends on the following parameter 
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where H(T = Mm) is the Hubble parameter when the temperature T becomes as low as 
the mass of the decaying right-handed neutrino, Mm- The so-called out-of-equilibrium 
condition is roughly given by Ki ^ 1, and the final lepton asymmetry is strongly suppressed 
for Ki ^> 1 [pO|l . The parameter Ki can be rewritten in terms of a mass parameter rnTj [^Ofi : 

ifj ~ mt , (35) 
O.OOleV V ; 

where 

We use the Yukawa couplings hik defined in the basis where the both mass matrices of 
right-handed neutrinos and light neutrinos are diagonal: 

(HI 



Y, hii hik Jf~. = mV i b ik ■ (37) 



Let us define a matrix Xij as follows: 



x v = I = ■ (38) 



Then from Eq. (|37|) one can show 

{x:<K) i . - 5, - {xx\ . (39) 
Thus, we see that the ffti parameter in Eq. ([36]) is bounded from below as follows: 

mi = J2 m k\ x ik\ 2 

k 

> mmim^YlXikl 2 



k 



Y X ik 

k 



> min{m^} 

= mm{m u A. (40) 

3 J 

In the last equation, we have used Eq. (|39|) . Therefore, if the light neutrinos are degenerate 
as m\ ~ m v 2 ~ ~ O(0.1)eV, the K parameter becomes at least as large as K > (9(100) 
[see Eq. fl35|) l, which leads to a strong suppression of the lepton asymmetry generated by 
the decays of thermally produced right-handed neutrinos. 
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Figure 1: Histograms of the ub/s, | U e 3 1 and Ai v . We randomly generate the small 
perturbations as in Eq. ([H]), and collect the data set of these quantities if the generated 
mass spectrum and mixing angles of the light neutrinos satisfy the conditions for the LMA 
solution in Eq. (|T4|) . Then we show the frequencies of these quantities in the vertical axes. 
The generated mass spectrum and mixing angles for the solar neutrino oscillation is also 
plotted in the lower-right figure, in which Amg 0lar is plotted with log 10 [Amg Olar /eV 2 ] unit. 
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Figure 2: The same as Fig.pl but for the LOW solution. 
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